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Exceptionally bright optical emission from a rare and dis- 
tant y—ray burst 
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Long y-ray bursts (GRBs) are produced by the dissipation of ultra-relativistic jets launched 
by newly-born black holes after the collapse of massive stars. Right after the luminous and 
highly variable y-ray emission, the multi-wavelength afterglow is released by the external 
dissipation of the jet in circumburst medium. We report the discovery of a very bright (~ 10 
mag) optical emission ~ 28 s after the explosion of the extremely luminous and energetic 
GRB 210619B located at redshift 1.937. Early multi-filter observations allowed us to witness 
the end of the shock wave propagation into the GRB ejecta. We observed the spectral tran- 
sition from a bright reverse to the forward shock emission, demonstrating that the early and 
late GRB multi-wavelength emission is originated from a very narrow jet propagating into 
an unusually rarefied interstellar medium. We also find evidence of an additional component 
of radiation, coming from the jet wings which is able explain the uncorrelated optical/X-ray 


emission. 
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On June 19, 2021, at 23:59:25 Universal Time (UT) the Swift Burst Alert Telescope trig- 
gered on GRB 210619B' and promptly distributed its coordinates. GRB 210619B was also de- 
tected by the Gravitational wave high energy Electromagnetic Counterpart All sky Monitor”, the 
Konus/WIND Experiment", and by the Fermi Gamma-Ray Burst Monitor”. This allowed ground- 
based robotic telescopes to rapidly repoint towards it, starting unfiltered optical observations while 
gamma-ray activity was still ongoing. D50^^ and FRAM-ORM' telescopes were able to capture 
the earliest optical emission. The data acquired since 74 + 28 s display an overall smooth optical 
light curve which faded from magnitude r — 10.6 to r — 17 about 3000 seconds later as shown in 
Fig. 1. Prompt observations by Mini-MegaTORTORA (MMT-9) camera” started at Ty 4- 55.5 s and 
provided high temporal resolution (1 s exposures) data and images of the transient in photometric 


Johnson B and V filters. 


These observations allowed us to probe the color evolution (see Fig. 2) of the GRB 210619B 
optical afterglow over the first hundreds of seconds since its onset. Following observations with 
the Gran Telescopio CANARIAS (GTC) revealed? the redshift of the burst, z — 1.937 placing the 
event among the most luminous (~ 2.5 x 10°erg/s) and most energetic ( 3.4 x 10% erg) ever 


detected (see Methods). 


The prompt emission which lasts for ~ 80 s is thought to originate from the internal dissi- 
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pation of the jet via shocks or magnetic re-connection in an optically thin region. The 


dominant contribution to the prompt emission is suggested to be the synchrotron radiation from 


non thermal population of electrons ! ^^ 


. To establish whether the observed early optical emis- 
sion is originated from the prompt emission, we performed time-resolved spectral analysis of the 
prompt emission using the time-bins corresponding to the data-sets provided by D50 and MMT-9. 
To model the GRB spectra in the widest possible energy range, we select the data provided by the 


Fermi/GBM instrument (8 keV - 10 MeV). For all the time-resolved spectra of Fermi/GBM we find 


that the optical flux exceeds the low-energy extrapolation of the GRB spectra by 2-2.5 orders of 
magnitude (see Fig. 1). This analysis excludes the possibility for the optical emission to arise from 
a single-component prompt emission spectrum ^. We also exclude that the optical emission could 
have been produced by an additional low-energy component in the prompt phase for the following 
reasons. In a standard optically thin emission regions, the bright optical emission would originate 
from a synchrotron cooling of non-thermal electrons, while the usual keV-MeV emission — as an 
Inverse Compton radiation from the same seed photons. In this model, the second order Inverse 
Compton component would be brighter than the keV-MeV emission by itself causing an energy 


budget problem '^"*. 


Alternatively, the 30-100 s optical emission could be the delayed prompt 
emission from the first bright episode of the GRB at 0-20 s. This would require a mechanism that 
can store some part of the energy of the prompt emission and release it later, at larger distances. 
High energy neutrons, primarily present in the fireball ^^^?! or produced on the prompt emission 
side, could in principle provide this mechanism by their later decay. However, the bulk Lorentz 
factor required in this scenario is rather too small I ~ 25” (for a delay of 50 s) compared to the 
Lorentz factor distribution of long GRBs”. 

A bright and primarily optical flash is expected from a reverse shock (RS, hereafter) propa- 
gating into the GRB ejecta ^^"^. The RS heats the particles of the ejecta and fades immediately 
after crossing them. Due to the energy sharing among densely populated particles, the average en- 
ergy gain of an electron is rather small compared to the forward shock scenario, making the peak 
of the RS emission to appear in the optical, infrared and radio bands. An order of — 20 optical 
counterparts of GRBs (see Fig. 4) and few late-time radio transients were suggested to originate 
from RS ^"^, Figure 4 shows the comparison among the early optical emission of GRB 210619B 
and all the bright optical emission interpreted as RS for other GRBs, along with the optical light 


curve of the highly variable naked-eye GRB 080319B '"^'. The RS candidate from GRB 210619B 


is the most luminous together with the historical GRB 990123". However, compared to GRB 


990123 which is characterised by few optical exposures, the optical light curve of GRB 210619B 


is the first observed in great details allowing us to accurately model it by external shocks. 


We first model the optical and X-ray light curves separately by series of power-laws (See 
Methods). The optical light curve is best represented by the RS component that peaks at ~ 30 
s with a following decay of t !9. Forward shock (FS, hereafter) is initially subdominant but it 
becomes later the only component in the optical band. The FS is best described by t^? before 
the break time of ~ 4 x 10% s and t !? later on. The X-ray light curve is simply composed 
by two power-laws smoothly connected at ~ 10* s with the shapes of t^"? and t~!° before and 
after. We preliminary interpret the optical/X-ray lightcurves as arising from the RS (only in the 
optical band) and the FS with an early jet-break at 10“ s. The temporal shapes and the optical/X- 
ray spectral indices (foy. ~ Dx © 0.9, where Fp optix X v ^) suggest that at least before the 
jet break the observed cooling frequency of the FS is above the X-ray band (0.5-10 keV). This is 
additionally supported by the measure of the Fermi/LAT (7100 MeV) spectral index of Spar = 
1.68 + 0.32. However, a detailed investigation of the shapes of the light curves suggests that 


the usual FS model***** is not capable to account for the chromatic X-ray and optical temporal 
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behaviour even if the jet structure or the presence of the second wider je is taken into 
account. To resolve the chromaticity problem, we invoke the Large Angle Emission (also called 
high latitude emission, LAE, hereafter) from a structured jet, i.e. delayed prompt like emission 


41,42 


from the jet envelope ^^, which accounts for an additional X-ray component (see Methods). 


The joint optical/X-ray light curves are further modelled by a combination of the RS, FS 
and LAE components (see Fig. 3). To account in details of forward shock model we use the 
publicly available code afterglowpy ^. We included empirically the RS component. The best fit 
model is described by the RS component peaked at ~ 36 s, with a temporal decay of x t "7 and 


a characteristic time of teu. Rs ~ 10? s, after which the RS component drops exponentially. The 


cutoff-time corresponds to the passage of the cooling frequency through the observed optical band. 
This allows us to constrain the RS cooling frequency ^ at the shell crossing time vV.rs(to,Rs) ~ 
3 x 1077 Hz. We clearly see the expected softening of the optical spectrum from few hundreds 
of seconds until ~ 10° s (Fig. 2). The best fit parameters of the FS component suggest that the 
cooling synchrotron frequency of the FS at the RS crossing time is ve rs(toRs) ~ 107% Hz. We 
then carefully estimate the relative magnetisation parameter*** Rg = 2.5 + 0.3 which gives the 
magnetic equipartition parameter of the RS component of eg rg = 0.12*920. In principle, the RS 
can be formed in a magnetised ejecta ^^"^. The FS model of GRB 210619B suggests that the 
circumburst medium is very rarefied, with a particle density of 6 x 10? cm ?. The densities as 


49,50,51,52,53 and 


low as 107? — 10^? cm”? were inferred from the modelling of a few GRB afterglows 
it could be connected with an unusual location of this GRB in the host galaxy or even be associated 


with an alternative progenitor such as a blue super-giant star””. 


The low circumburst medium implies that the jet was initially propagating almost in a vac- 
uum prior to the deceleration. Therefore, the ejecta has been spread and significantly decreased its 
magnetisation at the deceleration side providing the observed bright optical flash ^. This regime 
of formation of the RS emission (magnetised jet propagating into a vacuum ^) corresponds to 
egrs ~ Ty1/trs ~ 0.3 roughly consistent with eg rs = 0.124020, where we have assumed that 
the RS ejecta corresponds to the duration of the first bright prompt emission pulse (~ 10s). As 
predicted in this scenario, we observed equal deceleration times for the RS and FS. Since the mag- 
netisation of the RS ejecta has significantly decreased since times t < trs”, then the GRB jet 
could be initially highly magnetised, and can be even dominated by the Poynting flux ^^^". 

Rich optical/X-ray/GeV data allow us for the first time to fully constrain the contributions 


from the RS, FS and the LAE components. We have demonstrated that one needs to invoke a 


narrow (0.9*01 deg) structured jet propagating into a very rarefied medium (~ 6 x 107? cm~?) 


to explain multi-wavelength observations of the GRB 210619B afterglow. While this GRB is 
produced by rare (but realistic) combination of observables, it follows known spectral correla- 
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tions , and it shows a consistency with the standard internal shocks model for the production 


of the prompt emission (with an efficiency of 0.013+9:820) 663, Furthermore, its true kinetic energy 
Ekin.o = 3.824? x 10°? erg is consistent with the energy of the population of supernovae associ- 
ated with long GRBs. An opening angle of ~ 1 deg is rarely observed in GRBs” and is found 
in less than 10% of studied populations °° which are, however, usually embedded in environments 
with densities orders of magnitude larger. Extremely narrow opening angles < 1.5° and jets with 
steep structure k > 4 are expected in long GRBs due to the initial propagation of the jet into the 
progenitor material””””. However, the opening angles inferred from the observations of late-time 
afterglow steepening (jet breaks) suggest much wider jets^*. This discrepancy can be resolved by 
explaining the observed jet breaks as originated from viewing angle effects”. In this case, the 
"truly" on-axis GRBs are few and they expect to be the most luminous with the early jet break 
due to the narrow jet with a steep structure, exactly as witnessed in GRB 210619B. Importantly, 
the combination of observables suggests that the GRB jet could be initially magnetised”””””?. The 
fact that the GRB 210619B is located in a low-density medium and it is observed with the jet 
aligned towards us, allowed us to witness this rare and distant (z = 1.937) monster as source of an 


exceptionally bright optical transient (~ 10 mag). 
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Figure 1: Upper panel — early optical lightcurves from D50, FRAM-ORM and Mini- 
MegaTORTORA (channels with 1 s and 5 s temporal resolution) telescopes along with Swift-BAT 
count rate. Optical fluxes are corrected for extinction. Blue vertical bars denote the intervals for 
whose joint gamma-optical spectra are shown on lower panels. Lower panels — Fermi GBM spec- 
trum and its extrapolation towards optical range at different time intervals along with optical flux 
measured by D50. We also show the slope of the optical spectrum. The slope is estimated using 


data from several telescopes observing simultaneously in slightly different wavebands. 
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Figure 2: Upper panel — overall optical (D50) and X-ray (Swift XRT) light curves and their empir- 
ical decomposition into earlier reverse shock and later forward shock components. Middle panel — 
optical g — r color estimated from comparison of fluxes from simultaneous unfiltered observations 
by several telescopes (D50, FRAM-ORM, and Mini-MegaTORTORA, earlier points) and from 
direct measurements in photometric filters (D50 and data from GCN circulars, later points). The 
measurements are not corrected for Galactic and extragalactic extinction. Lower panel — corre- 
sponding F = vu”? spectral slopes after correcting for the Galactic extinction (with EG, — 0.17) 
and an additional extragalactic intervening dust absorber at z = 1.095 (with E, ,. = 0.40). 
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Figure 3: Optical (orange), X-ray (red) light curves together with the Fermi/LAT data points (ma- 
genta). We show the contributions from the reverse shock to the optical emission by dashed blue 
line and the forward shock emission by blue, red and magenta dotted line in the optical, X-ray and 
GeV bands, respectively. The Large Angle Emission is drawn by the dashed-dotted lines in the 
optical (blue) and X-ray (red) bands. The overall optical and X-ray models are shown by solid 
orange and red lines, respectively. 


Mean Posterior Values 


Parameters mode! 
RS+FS RS+FS+LAE 
Exin (7.3 40.7) x 1078 erg ani x 1098 erg 
6; 0.275, deg 0.9*93 deg 
k 1.044005 6.2 + 2.3 
2:005 ee x 10 ons. | DOBR x 1075 cm~’ 
Forward Shock x e el 
p 2031*55% 2.02350 6 
2 5.6+0:1 x 103 00575 
EB 0.381913 7.94474 x 1074 
trs 36 + 0.48 36 = 0.48 
Reverse Shock tRScut 1350: s 22901399 s 
ARS —1.74 0.01 —1.75 = 0.01 
R 7.9713 x 10% cm 
Large Angle Emission To 3607 209 
E, 2204319 keV 
Exin o (4.4 + 0.3) x 109 erg | 3.87177 x 109? erg 
E, 6 2.1105 x 10 erg | (4.4 1.4) x 1090 erg 
Inferred parameters Ny AT x 1073 13314107 
RB 20+1 2:5 £ 0.3 
€B,RS 67 0.12+9:20 


Table 1: Mean posterior values of the parameters of two tested models (with and without 
the Large Angle Emission). 
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Figure 4: Comparison between optical light curves of GRB 210619B (red solid line) and other 
GRBs with bright reverse shock component (except GRB 080319B). Two bright GRBs, namely, 
GRB 990123” and GRB 080319B "? are shown with blue and green solid lines, respectively. For 
GRB 080319B, data before and after 60 s are from the optical V-band and optical R-band, respec- 
tively. The grey solid lines represent the GRBs (GRB 021004 ”, GRB 021211 7',GRB 050525A *, 
GRB 061126”, GRB 081007 *, GRB 090102”, GRB 09102476777879, GRB 130427A*, GRB 
140102A*', GRB 140512A*, and GRB 190114C**) observed with (possible) RS component. 
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Materials and Methods 


Optical data 


D50° is a 50-cm telescope located at Ondfejov Observatory of Astronomical Institute of 
Czech Academy of Sciences, Czech Republic. It is equipped with Andor iXon Ultra 888 EMCCD 
camera and a set of photometric filters of the Sloan system. The telescope reacted to the trigger 
and started observing the position of the transient at 2021-06-19 23:59:53 UT (74 + 28 s, during 
the ongoing gamma-ray emission). In the following 2 hours, it acquired a series of unfiltered 10-s 
exposures, then followed by several longer exposures in the Sloan g, r, + and z filters at later times. 


1** data of an 


The data have been processed and calibrated photometrically using PanSTARRS DR 
ensemble of field stars, simultaneously deriving both the zero point and color term for unfiltered 
images. The best fit for the instrumental photometric system is mag = r — 0.05(g — r) + ZP, very 


close to PanSTARRS r. 


FRAM-ORM' is a 25 cm f/6.3 telescope located at Observatorio del Roque de los Mucha- 
chos, La Palma, Spain. The telescope is equipped with B, V, R and z filters, and a custom Moravian 
Instruments G2-1000BI camera based on a back-illuminated CCD47-40 chip. It started observing 
the transient position at 7,+30 s and acquired a series of unfiltered 20-s exposures until 704-1050 s. 
Due to strong wind at the telescope site the images were distorted. In order to exclude the possible 
varying influence of nearby stars on the photometric measurements, the frames have been passed 
through a difference imaging pipeline based on HOTPANTS image subtraction code ^. The mea- 
surements have been calibrated using field stars and PanSTARRS DRI catalogue". The best fit 
for the instrumental photometric system is mag = r +0.10(g — r) + ZP, very close to PanSTARRS 


r, although slightly different from the system of D50. 


Mini-MegaTORTORA (MMT-9) nine-channel wide-field monitoring system" is a set of 9 
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7-cm objectives equipped with Andor Neo sCMOS cameras and independently installable Johnson- 
Cousins B, V, R and polarimetric filters, and mounted on 5 mounts. The system responded to the 
BAT trigger and observed the position of the transient since 2021-06-20 00:00:20 UT (70455 s) and 
until 2021-06-20 00:10:28 UT (7, + 663 s). The system simultaneously acquired series of frames 
with 1 s exposures, 5 s exposures and 30 s exposures in white light, 10 s exposures in B filter, and 
10 s exposures in V filter. The transient is clearly detectable in all the acquired sequences except in 
the B filter ones. As the field was crowded due to the large pixel scale of Mini-MegaTORTORA, 
the photometry has been performed on the images after the subtraction of a template image us- 
ing the HOTPANTS image subtraction code ^. The template has been constructed by averaging a 
large number of frames from observations of the same position on successive nights. Photometric 
measurements in Johnson B and V filters have been calibrated to Vega magnitudes using a large 
ensemble of field stars and their photometric measurements in Gaia DR2 catalogue converted to 
Johnson-Cousins system using conversion formulae based on Stetson standards. Also, both white 
light images and the ones in V filters have been calibrated to AB magnitudes using PanSTARRS 
DRI catalogue ** data of field stars by fitting the color terms and zero points for every sequence. 
The best fit for the unfiltered photometric system is mag = r + 0.353(g — r) + ZP, while for the 
V filter is mag = r + 0.538(g — r) + ZP. 


High temporal resolution Mini-MegaTORTORA data are acquired from the channels 1 and 
2 that were operating with 1 and 5 s exposures. Their light curves are generally similar to the 
D50 one acquired with 10 s temporal resolution, and do not show any significant excess variability 
relative to its smooth behaviour (see Fig. S1), despite still ongoing and highly irregular gamma-ray 
activity. This excludes that optical and prompt gamma-ray emissions are produced by the same 


process, and it is consistent with optical emission being mostly produced by a reverse shock. 


Early optical color has been estimated using simultaneous observations in slightly different 
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photometric systems. If both systems are represented as r-- C(g —r), then the ratio of fluxes among 
them is F; / F, = 10~°-4*(C1-C2)(9-") | We computed the average ratios of fluxes from all early data 
sets (unfiltered FRAM-ORM, unfiltered Mini-MegaTORTORA, V filter Mini-MegaTORTORA) 
relative to the interpolated flux from D50 in several temporal intervals. We then fitted them to 


evaluate the color index. The results are shown in Figure 2 and Table S3. 


Galactic extinction towards the transient has been taken from NASA/IPAC Extragalactic 
Database (NED) Extinction Calculator, and then used to correct the optical fluxes used for plotting 
Fig. 2. It corresponds to Ay = 0.467 and E? „ = 0.17. Later time spectral’ and multicolor** 
data suggest the presence of an additional reddening of the afterglow, shown as a strong broad 
absorption feature at ~4500 AA supposedly corresponding to a 2175 AA dust feature due to an 
intervening absorber at z = 1.095. To model it, we used the data" in g, r, i and z filters measured 
at around 74 + 10000 s, and we assumed a power-law optical spectrum and a Milky Way like 
absorption law"? with Ay = 0.30 + 0.05. This results in E; , = 0.40. The fit required no 


additional reddening due to the burst host galaxy per se, so we assume that it is negligible. 


X-ray data We have retrieved the X-ray (0.3-10 keV) lightcurve of GRB 210619B from the Swift 
Science Data Center supported by the University of Leicester". Four time-bins in the Window 
Timing mode (335-1864 s) and ten time-bins in the Photon Counting mode (209-283 s ,6395- 
1.3 x 10° s) were selected for the spectral analysis. The source and background spectral files, 
the ancillary response and the response matrix files were produced by the Online Swift-XRT GRB 
spectrum repository ". We have fitted all 14 unbinned XRT spectra by a simple power-law model 
using XSPEC (12.10.1) and adopting C statistic likelihood. To account for the Galactic”! and 
intrinsic metal absoption we applied tbabs and ztbabs models, respectively. The hydrogen column 
of the host galaxy Ny has been left as a free parameter. We do not find significant variations of Ny 


among 14 spectral fits. Finally, we have estimated unabsorbed X-ray flux in the 0.5-10 keV band 
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and used it for modelling the joint X-ray/optical afterglow lightcurves. We report the flux of all 14 


time-bin in the Table S1 in Supplementary Materials. 


Fermi/GBM data We have extracted the time-resolved spectra of Fermi/GBM corresponding to 
the time-windows of D50 and MMT-9 observation times. Fermi/GBM is composed by 12 sodium 
iodide (Nal) and two bismuth germanate (BGO) scintillation detectors". For each time-bin we 
have chosen spectral data from two Nal detectors and one from a BGO detector. The choice of 
Nal and BGO detectors is based on the source position angle. We have chosen the Nal-4 (51 deg), 
Nal-8 (26 deg) and BGO-1 (76 deg) detectors. We have extracted the Fermi/GBM spectra by the 
GTBURST tool. To model the background, we have fitted pre- and post-burst data by the energy- 
and time-dependent polynomial. The time-resolved spectra have been fitted by Band function ^? 
using XSPEC (12.10.1) and applying the PGSTAT likelihood. We took the spectral data of Nal 
from 8 keV to 900 keV, and the BGO data from 300 keV to 10 MeV. The best fit parameters are 
listed in the Table S2 in Supplementary Materials. We report the low energy photon index a, the 
high energy photon index 9, the break energy Ep and the observed flux in the range between 8 
keV and 10 MeV. In order to estimate the isotropic equivalent energy Eiso and the peak luminosity 
Lp iso Of GRB 210619B, we also model with the Band function the time-averaged spectrum and 
one-second peak spectrum. The time-integrated spectrum is defined as of the duration of the GRB 
(Togo = 54.785 s from 0.576 s to 55.361 s as defined by the Fermi/GBM Burst catalog). The one- 
second peak time spectrum is defined between 0.51 and 1.51 s as reported by Fermi/GBM team’. 
Due to the extreme brightness of the spectrum, we notice the presence of the Iodine K-edge at 
33.17 keV"? and therefore we ignore the Nal data from 30 keV to 40 keV for the time-integrated 


and one-second peak time spectrum. We also extend the spectral coverage by BGO up to 40 MeV. 


GRB 210619B and spectral-energy correlations The best fit parameters of the time-integrated 


spectrum are the following: the low-energy photon index a = —0.86 + 0.02, the high-energy 
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photon index 8 = —2.1879:93, the break energy in the photon spectrum Eg = 178.67* 715 keV 
with PGSTAT equal to 788 for 339 degrees of freedom. The flux of Too spectrum is estimated 
to be (6.67 + 0.08) x 10~°erg cm~*s~! between 0.34 keV and 3404.8 MeV. This band is used 
to compute isotropic equivalent energy Eiso and the isotropic equivalent peak luminosity Ly ;so 
in the energy band 1 keV - 10 MeV in the rest frame given the redshift of 1.937. We estimate 
Eiso = (3.44 + 0.04) x 10% erg, and the peak energy of the energy flux spectrum in the rest frame 
Epeak = 5987133 keV. The total energy and the spectral peak are in agreement with the Amati 


relation "(see Fig. S3). GRB 210619B is the second most energetic GRB after GRB 160625B. 


The best fit parameters of the one-second peak time spectrum are the following: the low- 
energy photon index a = —0.36 + 0.03, the high-energy photon index 5 = —2.16 + 0.03, the 
break energy in the photon spectrum Ey = 212.88* 112; keV with PGSTAT equal to 548 for 339 
degrees of freedom. We estimate the flux of one second spectrum (8.8 + 0.1) x 10 ?ergem ?s^! 
between 0.34 keV and 3404.8 MeV. The corresponding isotropic equivalent peak luminosity in the 
rest frame is Ly iso = (2.45 £0.03) x 10% erg s^! and the peak energy of the energy flux spectrum 
in the rest frame is Epeak = 1024*77 keV. The peak luminosity and the corresponding peak energy 
are consistent with the Yonetoku relation "(see Fig. S4). GRB 210619B is positioned as the second 
most luminous GRB after GRB 110918A and it has the largest spectral peak energy among long 


GRBs in the Yonetoku plane. 


Fermi/LAT 


The Large area Telescope (LAT) on board Fermi consists of a tracker and a calorimeter sensi- 
tive to the gamma-ray photons of energy between 30 MeV and 300 GeV '**. We used GTBURST 
tool from the official Fermi-software to extract and analyse the data. Due to its position (R.A. 


= 319.71? and Dec. = 33.86"), the source was out of the field of view (FoV) of LAT during 
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the first 200 seconds after trigger. The LAT data were analyzed in two periods: a) To + 200s 
to To + 1800s, and b) 7, + 6400s to 7, + 10000s in the energy band of above 100 MeV. For 
this analysis, we used a region of interest (ROI) of 12° centred at the burst position provided by 
Swift/XRT ^. In addition, a zenith angle cut of 100° was applied to reduce the contamination of 
the gamma-ray photons from the Earth limb. We assumed a spectral model of type "powerlaw2" 
for the source. We considered "isotr template" and "template (fixed norm.)” for the particle back- 
ground and the Galactic component, respectively. The calculations of the flux in the energy bin 
given above are calculated with the "unbinned likelihood analysis" considering a minimum test 


statistics (TSmin) of 10. The Fermi/LAT analysis shows that for period a) the source is found with 


an energy flux of (2.843-2.70) x 10? erg s^! cm ? and with a TS value of 12. The spectral index 
for this period is -1.68+0.32. For period b), the TS of the source is 8 and the flux upper limit is 


9.13x10 ^7? erg s! cm ?. 


The empirical model for the X-ray/optical light curves We first model the X-ray (0.5-10 keV) 
and optical lightcurves empirically in order to assess the temporal evolution of the possible reverse 
and forward shock components. The X-ray lightcurve is modelled by a smoothly connected broken 


power-law function 


—(1/n) 
1 t "BRP M t Bii 
Ka) SPORE | a | 
v Fux(t) = Fo 2 \ toreak,FS,X 2 \ toreak,FS,X u 


which returns an initial decay of the energy flux vF* œ t~°° before the temporal break of 


break Fs,x ~ 10^ s. This is then followed by a flux decay vF,x x t 1? 


after the temporal break. 
We have adopted an intermediate smoothness parameter n — 5. The spectral index in 0.5-10 keV 
remains roughly the same £x = 0.9 (where F, x « v—°x) before and after the temporal break (see 


Fig. S6 for the distribution of the parameters). 


To account for the complex temporal structure of the optical lightcurve, we adopt two smoothly 


broken power-law functions, one representing the reverse shock component (RS hereafter) and the 
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other the forward shock component (FS hereafter). Since we do not observe the rise of the early 
optical emission, we fix the initial rise index of the RS component to the theoretically driven value 
of airs = 0.5?*, where FL Rs,opt X t^^^5, To account for the adiabatic cooling of the ejecta ^ 
which has previously experienced the passage of the RS, we introduce an exponential cutoff in the 
empirical function of the RS at t, 5s. The candidate FS component is modelled as in the case of 
X-ray lightcurve, with the only difference that we use F, opt data as input for the empirical function 
1. The model returns a peak flux density of the RS component of F985 ~ 0.5Jy and a decay index 
of aa rs = —1.6, a pre-break decay index of the candidate FS component of o psy. ~ —0.6, a 
temporal break of tprcak,FS,opt Z 4 X 10‘s and a post-break index of Q2.FS,opt © —1.2. The cut- 
off time teut rs Of the RS component is not constrained in the fitting by the empirical model (see 


Fig. S7 for the distribution of the parameters). 


The spectral index in the R-band is derived by the color differences throughout all the ob- 
served emission. It starts (at ~ 30s) with a quite low value of Bope = —2.10 (however with huge 
uncertainty of 2.43), then it fluctuates around bopt = 1.0 rising towards softer spectral values until 
the end of the RS component. With the rise of the candidate FS component the optical spectral 
index shows a softening trend from a value of Bope = 0.37 + 0.26 to 5,,, = 1.40 + 0.32. Since the 
redenning from the host galaxy is unknown, the absolute values £,,, is uncertain but the overall 
softening of the spectra during the RS and FS components is expected to remain the same, since 
the additional correction by the host would affect all the values by the same ô bopt. It is remarkable 
the overall consistency of the spectral indices derived accurately in the X-ray band and by the color 
differences in the optical bands, except that of ~ 1000s, at the time of transition between RS and 


the candidate FS components. 


We then compare the spectral shapes and the temporal behaviour of optical and X-ray emis- 


sion with the expectations from the standard afterglow model of GRBs****. The standard FS 


24 


model is based on the ultra-relativistic, adiabatic, self-similar spherical solution of Blandford & 


McKee for a dynamics of a point like explosion”. 


The radiation is computed assuming a syn- 
chrotron emission from shocked-accelerated non-thermal electrons. For our analytical estimates, 
we ignore the effects of the synchrotron self Compton (SSC) cooling of electrons. This system is 
characterised by the following usual parameters: Ex, the kinetic energy of the spherical ejecta, n 
the density of the circumburst medium, e, the fraction of kinetic energy given by non-thermal elec- 
trons, eg the fraction of the kinetic energy in the magnetic field, p the spectral index of power-law 
distributed electrons, and 6; the opening angle of the jet. We assume that all electrons that received 
€, Eja, got accelerated. We also assume an homogeneous circumburst medium. We further argue 
in terms of the synchrotron Vm and cooling ve frequencies, which are the observed frequencies cor- 
responding to the minimum Lorentz factor of accelerated electrons and the cooling Lorentz factor 


of electrons, respectively. We further compare the observed temporal and spectral indices with 


expectations from the FS model”®”’, 


We assume that before the tpreak,rs,x the X-ray band is located between Vm and ve. Consid- 
ering an homogeneous circumburst medium, in order to produce the initial decay of vF* œ t^ 99, 
one would require an electron index of p ~ 2.2 which would imply an XRT spectral index of 
Bx = 0.6 harder than observed £x oy, © 0.9. In our assumed case, i.e. when Vm < vx < v, at 
t < tpreak,rs,x, the temporal break tpreak,Fs,x could be caused by the passage of 1, through the XRT 
energy band at tprcak,Fs,K- However, it is highly inconsistent with the observed temporal break of 
the FS component in the optical light curve tpreak,FS,opt © 4 X 10? s, since the expected time of this 
spectral transition in the R band is tpreak,v.,opt > break ps x (75 EE )2 > 2 x 10" s. Therefore, 
we can exclude the scenario when the temporal optical/X-ray break is associated with the cooling 


frequency. 


We are then left with another interpretation of the temporal break, which is an early jet 
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break ^*^? at ~ 10s. To identify the temporal shape of the light curve after tprcak,FS,K © toreak FS opt 
we first assume that Vm < Vopt,x < Ve at ~ 10?s. This is justified by the rough consistency of 
the observed optical and X-ray spectral indices Bx 7 opt- For the previously inferred p = 2.2, 
the temporal index after the ~ 104s should be ax = Qopt = 3p/4 ~ 1.65 steeper than observed 
VF,x x t "9. However, if we assume that Vm < vg < Ve while vx > ve, we would require very 
soft p = 1.9. And if we insist on p = 2, the temporal index after the jet break would be even 
steeper ay = 1.75. Therefore, the most optimal scenario would correspond to the case of vx S ve, 
i.e. the cooling frequency is slightly above the XRT energy range. This is additionally supported 
by the quite hard value of the spectral index measured in the Fermi/LAT band (30 MeV-300 GeV) 
Butar = 0.68 + 0.32 which corresponds either to the very end of the synchrotron spectrum or to the 
initial rise of the SSC component. While the most optimal scenario (vx € ve) is roughly consistent 
with the temporal shape of the X-ray light curve, it is highly inconsistent with the joint chromatic 
optical light curve. The optical light curve is shallower than the X-ray light curve before the tem- 
poral break (01 ps opt © —0.6 VS Qi psx 7 —0.9) and after (o Fg opt © —1.2 VS as Rg x S —1.5). 
Additionally, the observed temporal breaks are not simultaneous as expected from the jet break and 
the X-ray spectra are much softer than expected in Vm < Vopt,x < Ve. Therefore, one needs to in- 
voke an additional process/component in order to describe both optical and X-ray light curves. We 
also notice that the wind medium would return more problematic solutions for the FS component, 


since even for Vm < Vopt,x < Ve, very soft electron spectrum p = 1.5 would be required. 


X-ray/optical light curves from external shocks We model the joint optical and the X-ray light 
curves with the FS model using publicly available code afterglowpy ^. The adopted model corre- 
sponds to a Power law structured jet (without the spreading effects) decelerating in a homogeneous 
circumburst medium. We have selected a power-law structured jet because we expect to see the dif- 
ference in the shape of the light curve compared to the top-hat jet at late times if the early temporal 


break at ~ 10^ s is associated with the jet break. We set the viewing angle of the observer to zero, 
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since we observe one of the most luminous and energetic GRB. We include the synchrotron and 
SSC radiation components. To account for the RS component, we adopt the previously described 
phenomenological function with fixed rising shape and the smoothness parameter. We exclude the 
fractional parameter x (fraction of accelerated parameters) since we found that it is not constrained 


and it does not change the typical values of other parameters. 


The modelling returns a jet with an isotropic equivalent of kinetic energy of Exi, ~ 7 x 1056 
erg and the opening angle of ~ 0.2 deg, propagating in a rarefied circumburst medium with n ~ 
2 x 10 ?cm ?. The microphysical parameters are the following: e. ~ 0.006, eg ~ 0.4 and 
p = 2.03. In this joint X-ray/optical modelling we were able to better constrain the parameters 
of the RS component: tors ~ 36s, FSRS ~ 0.5 Jy, aaps © —1.7, teurs S 1.4 x 10?s. These 
parameters are obtained by considering the optical data corrected for the Galactic absorption and 
the absorber at z=1.1, but it does not take into account the unknown reddening by the host galaxy. 
The combined FS and the RS models together with the optical/X-ray data are shown in Fig. S2. 
One can notice that the FS model is not able to reproduce the observed chromatic light curves, 
as expected from the previous analytical estimates. We have additionally tried to model the light 
curves by the power-law structured jet, two-component jets and we have failed to reproduce the 
fast declining X-ray light curve together with chromatic shallower optical light curve. One could 
assume an additional energy injection to the decelerating forward shock '?”!!! to flatten the optical 
light curve, however this would necessarily produce similar temporal behaviour also in X-rays '? 


which is not observed. 


X-ray/optical light curves from external shocks and the Large Angle Emission To account 
for the chromatic X-ray and optical temporal behaviour, we have included the Large Angle Emis- 
sion (LAE)'?”!?? i.e. the delayed prompt like emission received from the higher latitudes of the 


jet. It was suggested that the LAE from structured jets could account for the diverse X-ray after- 
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glows observed in GRBs, including plateaus, normal decays (FS like) and late-time exponential 
declines ^^. We have considered LAE from the jet wings, i.e. the regions outside of the jet core. 
The structure of the bulk Lorentz factor and the angular energy distribution follows a power-law 
profile with index k, i.e. ['(@),€(0) x 87* for 0 > 0;. We have fixed the comoving spectrum to 
the Band function with a = —1.0 and 6 = —2.50. Thus, the parameters that define the LAE are 
the following: the bulk Lorentz factor of the jet head Io, the transparency radius R (the size of the 
jet head where jet wings are optically thin ^), E, the "observed" peak energy of the emission, the 
jet opening angle 6; and the jet structure index k. Under our assumptions, the jet opening angle 6; 
and the power-law profile index k are common parameters for the LAE and FS models. We have 
included the FS and RS models together with the LAE self-consistently. The combined RS, FS and 
LAE model is capable of describing the complex optical and X-ray light curves. The early X-ray 
lightcurve is initially dominated by the LAE, while the optical light curve is entirely made by the 
early RS and later on FS components only. The jet of GRB 210619B is best described by the fol- 
lowing parameters: the isotropic equivalent of the kinetic energy Epin = 2.5* 19! x 1096 erg, the jet 
opening angle of 0; = 0.9*5 deg, the power-law index of the structure k = 6 + 2. It corresponds 
to the opening angle-corrected prompt emission energy of Eg = 4.4 + 1.4 x 10°° erg (consistent 
with the Ghirlanda relation®'), the true kinetic energy of the jet of Ej, = 3.8134? x 105? erg and 
the efficiency of the prompt emission of 0.013*0 029. The micro-physical parameters of the FS are 
p = 2.02370 008. eers = 0.051303 and egrs = 7.917" x 1074. The density of the circumburst 
medium is n = 5.5*19 x 107? em ?. The RS component is peaked at 36 + 0.4 s and declines 
with a temporal index of —1.75 + 0.01 with the following exponential cutoff at (2.3*03) x 10? s. 
The temporal decline of the RS?" suggests the electron distribution index of prs = 2.0. The LAE 
component that fits most of the X-ray emission is best represented by the following parameters: 
the transparency radius of the jet core R = (8*7) x 10? cm, the jet core bulk Lorentz factor of 
To = 360789) and the peak energy of the spectrum in the observed frame highly unconstrained and 


is spanning over the range of ~ 10 — 5600 keV. 


28 


The magnetisation of the jet The ratio between the cooling freguencies of the RS and FS compo- 
nents at the shell crossing time, i.e. at to,rs, defines roughly the relative magnetisation parameter 
Ve RS/Vc FS © Rg t where Rg = €p,ns/ egrs” . The cutoff time suggests that the cooling fre- 


quency crosses the observed R-band at ~ 2.3 x 10° s. We estimate the cooling frequency of the RS 


to,RS 
Ücut,RS 


—54/35 
) ~ 3 x 10!" Hz (conservatively, in the 


component at the shell crossing time ^ as vg ( 
thin shell scenario). Given the best parameters for the FS, we estimate the the cooling freguency at 
to,Rs AS Ve ps(to,RS) ~ 10/5 Hz. Therefore, the relative magnetisation parameter is Rp = 2.5 + 0.3 


which returns the magnetic equipartition parameter of the RS ep ps = 0.12+9:20, 


Markov Chain Monte Carlo sampling We have used the Goodman & Weare’s Affine Invariant 
Markov chain Monte Carlo (MCMC) Ensemble sampler ^ implemented by the python package 
emcee!” to explore the parameters of our empirical and physical models. We use the sums of 
Gaussian log-likelihoods in order to self-consistently estimate the model parameters when consid- 
ering optical and X-ray data simultaneously. We allow the MCMC to run until the number of steps 
exceeds fifty times of the maximum of the auto-correlation time. We report the 16th, 50th, and 


84th percentiles of the samples in the marginalized distributions throughout this paper. 
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Figure S1: Relative flux variations of the high temporal resolution (1 s and 5 s exposures) data 
from Mini-MegaTORTORA. The flux is normalized to the smooth curve corresponding to D50 
data. 
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Figure S2: The FS+RS model and the optical/X-ray/LAT data. 
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Figure S3: Amati relation for GRB samples with known redshift and estimated Epeak. The quan- 
tities Ej, and Eiso represent the rest frame peak energy and the isotropic energy, respectively. The 
shaded region marks the 3c scatter of the point distribution around the best fit line'?”. The blue 
markers represent the 46 GRBs detected with known E, and redshift and retrieved from '?" and 
references therein. GRB 210619B (this work) belongs to the 30 contour around the Amati corre- 
lation, and is one of the brightest sources. GRB 130427A !™ is marked with a green square which 
showed largest isotropic energy releases. The historical outlier, GRB 160625B ''', is shown with a 
red square. GRB 110918A '?? with the highest isotropic-equivalent luminosity observed so far, is 
marked with a purple square marker. 
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Figure S4: Yonetoku relation, where quantities L, and E, represent the 1-second peak luminosity 
and the rest frame peak energy, respectively. The data set '* (blue markers) is a collection of short 
GRBs, long GRBs and GRBs with high redshift (z>6). The shaded region marks the 3o scatter 
of the point distribution around the best fit line". GRB 210619B (this work) belongs to the 30 
contour around the Yonetoku correlation. The GRBs, namely GRB 130427A "8, GRB 160625B '”’, 
GRB 110918A '” are marked with green, red and purple square markers, respectively. 


35 


Qu 


Q92 


M MI 


logiotoreak [s] O4 Wo 


Figure S5: The joint corner plot for the parameters of the empirical model for the X-ray/optical 
light curves. 
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Figure S6: The corner plot for the parameters of the empirical model for the X-ray light curve. 
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Figure S7: The corner plot for the parameters of the empirical model for the optical light curve. 
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Figure S8: The corner plot for the parameters of the model which includes the emission from the 
reverse shock and the forward shock only. 
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Figure S9: The corner plot for the parameters of the full model including the emission from the 
reverse shock, the forward shock and the Large Angle Emission. 
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time-bin observing mode Flux 


(s) (10-72 erg cm~? s^!) 
209-283 PC 1504.771730] 
335-500 WT 897.727 9159 
500-1000 WT 555.6640183 
1000-1500 WT 32306050 
1500-1864 WT 283,84! 1162 
6395-7565 PC 58.907502 
11828-13237 PC BAT. 
17763-19024 PC DAMITO 8 
23610-34828 PC 8.08*0011 
50936-58400 PC 2002 06 
58400-99305 PC 22076 
109476-196624 PC [045505 
213000-400000 PC 0260 


Table S1: Unabsorbed XRT Flux [0.5-10 keV] 


time-bin Optical instrument a B Eo Flux 10% — PGSTAT/d.o.f. 
S keV erg cm~? s7! 


55.5-56.5 MMT-9 [channel 1] —0.86+9:23 —2.667973 63.572906  2.00*433 391/318 
58.7-59.7 MMT-9 [channel 1] —0.38'072; —2.64+915 34.332057  2.50+9:25 352/318 


1 
60.9-61.9 MMT-9 [channel i] —0.83*57; -—10.00*900 45.841325 1.447003 382/318 


NEUES NETS ETE 392/318 
EE FI SIE un 399/318 


55.5-60.5 MMT-9 [channel 2] —0.8710-° —2.59*01 
60.6-65.6 MMT-9 [channel 2] —0.69*925 —2.69*$: 


28.2-38.2 D50 —0.95t997 280. 13660" 241072; 370/318 

39.3-49.3 D50 —1.0119% 2206 119,554+1234 93 97+0-18 473/318 
0.05 0.20 11.09 0.17 

50.8-60.8 D50 —0.91+9:07 —2.43+0:09 18h20 p 8.082540 522/318 

62.0-72.0 D50 —0.75+039 2024373 9455713840. (491 9:09 388/318 
0.67 0.30 7.54 0.05 


Table S2: The best fit parameters of the time-resolved Fermi/GBM spectra modelled by 
the Band function. The flux is computed on the range from 8 keV to 10 MeV. 
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Interval FRAM Clear MMT-9 Clear MMT-9V  MMT-9B (g — r) B 
S F/ Fpso Ff Fpso Ff Fpso Ff Fpso 
30 .. 50 1.01 £ 0.1 —0.06 + 0.75 —2.10+ 2.43 
50 .. 100 0.85 + 0.01 0.71 +0.03 0.61 + 0.05 < 0.18 0.99 + 0.07 1.42 + 0.26 
100 .. 200 0.88 + 0.01 0.76 +0.06 0.62 €: 0.10 < 0.13 0.90 + 0.09 1.09 + 0.33 
200 .. 400 0.88 + 0.01 0.80 +0.12 0.55 Æ 0.20 < 0.36 0.88 + 0.10 1.08 + 0.37 
400 .. 600 0.90 + 0.02 0.15+0.12 0.96+0.44 < 0.65 0.99 + 0.14 1.53 £ 0.46 
600 .. 1000 0.83 = 0.02 1.34 £ 0.19 2.56 £ 0.62 
Late-time data 

3200 0.67 + 0.07 0.37 + 0.26 
5500 0.57 £ 0.13 0.00 + 0.47 
8640 0.76 + 0.04 0.70 + 0.16 
13824 0.82 + 0.04 0.92 + 0.16 
66528 0.95 + 0.09 1.40 + 0.32 


Table S3: Average fractions of fluxes from FRAM and MMT-9 relative to D50 in different 
time intervals, and corresponding (g — r) color indices and spectral slopes 5. The fluxes 
and color indices are not corrected for the Galactic extinction, which contributes EC , = 
0.17 in the direction of the burst. The 8 spectral slopes, on the other hand, are corrected for 
both Galactic extinction and an additional E, ,. = 0.40 dust extinction due to an intervening 
absorber at z = 1.095. For completeness, the colors and the corresponding slopes derived 
from the data in Tables S4 and S5 are shown. 
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Time, s Mag Exposure, s Filter 


2771 17.80 0.05 g 131 
2902 17.93 + 0.06 g 123 
3061 17.87 £0.05 g 188 
3442 17.20 + 0.05 r 60 
3505 17.20 + 0.06 r 60 
3569 17.19+ 0.06 r 60 
3633 17.30 + 0.06 r 60 
3696 17.30 + 0.07 r 60 
3760 17.29 + 0.06 r 60 
3823 17.23 + 0.06 r 60 
3886 17.27 + 0.06 r 60 
3952 16.88 + 0.05 a 60 
4014 . 17.05 + 0.07 a 60 
4078 16.93 + 0.06 a 60 
4142 . 17.00 + 0.07 a 60 
4206 16.91 + 0.07 a 60 
4271 | 17.00 + 0.07 a 60 
4335 . 17.03 + 0.07 a 60 
4399 17.01 +£ 0.08 a 60 
4462 . 17.02 + 0.08 a 60 
4589 16.77 + 0.06 z 186 
4781 16.81 0.08 Ź 188 
4971 16.68 + 0.06 z 188 
5384 18.14 0.10 g 632 
5799 17.57 £ 0.08 r 186 
5988 17.57 0.09 r 186 
6177 17.56 0.10 r 187 
6865 16.90 £ 0.11 z 1178 
84656 20.16 + 0.19 r 180 
163842 21.12 + 0.48 r 180 


Table S4: D50 data in photometric filters 
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Time Mag Filter Reference 
9498 17.90-50.20 V  GCN30265% 
21492 18.51+0.10 R GCN 30268?” 
19824  18.51+0.10 R  GCN30268* 
9490 17.44 + 0.04 i GCN 30271 °° 
13848  18.11+0.03 r GCN 302715? 
9307 18.58-50.00 g GCN 302712? 
9399  17.82+0.03 r GCN 302712? 
14149 17474004 z GCN 302715? 
13939  17.74+0.04 i GCN 302712? 
9700  17.17+0.04 z GCN 302712? 
13755  18.93+0.03 g GCN 302715? 
39852  18.80+0.10 N  GCN30273'% 
28980  18.50+0.10 N  GCN30273'% 
51948 19.00+0.10 R  GCN30277'° 
13188 16.10+0.10 r  GCN30280''° 
72180 21.14+0.07 g | GCN30286'" 
65700 21.13+0.07 g  GCN30286'"' 
67284  20.18+0.05 r  GCN30286''' 
68904 20.09+0.06 4  GCN30286'' 
22104  18.52+0.12 r  GCN30288'* 
85679 19.80-20.10 R  GCN30291'? 
84154 20.50+0.40 r  GCN30292''4 
81215  19.90+0.20 R  GCN30293''° 
103699 18.88+0.11 z | GCN30294''® 
104881 19.65+0.02 i GCN 30294 ''* 
102244 19.92+0.02 R  GCN30294''° 
171876 20.60+0.30 R  GCN30299''7 
171289 20.730.008 R  GCN30303''® 
99792 20.000+0.20 R  GCN30305''9 
255989 21.40+0.10 R | GCN30309 ^ 
256068 20.09+0.28 i GCN 30320 1%! 
256068 21.45+0.33 r | GCN30320 ^ 
529067 22.56+0.16 r | GCN30338 ^ 
69088 19.47+0.06 R  GCN30791'* 
159856 20.43+0.05 R  GCN 30791? 
236914 21.82+0.26 r  GCN30791'*° 
241580 21.21+0.09 R  GCN30791 ^ 
593153 22.554+0.24 R | GCN30791 ?? 
1111648 23.50+0.30 R  GCN30791 '*° 


Table S5: Photometric data collected from GCN Circulars 
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